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Abstract 

For standard leptogenesis from the decay of singlet right-handed neutrinos, we 
derive source terms for the lepton asymmetry that are present in a finite density 
background but absent in the vacuum. These arise from cuts through the vertex 
correction to the decay asymmetry, where in the loop either the Higgs boson and 
the right-handed neutrino or the left-handed lepton and the right-handed neutrino 
are simultaneously on shell. We evaluate the source terms numerically and use 
them to calculate the lepton asymmetry for illustrative points in parameter space, 
where we consider only two right-handed neutrinos for simplicity. Compared to 
calculations where only the standard cut through the propagators of left-handed 
lepton and Higgs boson is included, sizable corrections arise when the masses of 
the right-handed neutrinos are of the same order, but the new sources are found 
to be most relevant when the decaying right-handed neutrino is heavier than the 
one in the loop. In that situation, they can yield the dominant contribution to the 
lepton asymmetry. 



1 Introduction 



Leptogenesis [l[ is often studied in parametric regimes where the masses of the right- 
handed neutrinos are either hierarchical or degenerate. The hierarchical limit is par- 
ticularly useful for gaining valuable analytical insights into the connections between 
leptogenesis and the observed neutrino oscillations j2|-Q]. To be specific, we discuss 
here the simple case of two right-handed neutrinos N\ }2 with masses M 12 and Yukawa 
couplings Y 1>2 to the Higgs and lepton doublets, i.e. a model as specified in Ref. 
The key simplification in the hierarchical limit, M\ <C M 2 , is that the evolution of 
the lepton asymmetry during leptogenesis as a function of z = M\jT depends up to a 
proportionality factor only on the ratio Mi/(Y 1 2 mpi), which characterises the washout 
strength. Here, T is the temperature and mpi is the Planck mass. The remaining pro- 
portionality factor characterising the amount of CP violation is Im[Y" 1 2 Y 2 * 2 ]M 1 /M 2 , up to 
corrections of order Mf/M|, which one neglects in the hierarchical limit. On the other 
hand, mass-degenerate right-handed neutrinos lead to resonant leptogenesis [6-lQ|. This 
corresponds to a phenomenologically attractive scenario, since the decay asymmetry is 
enhanced, which allows for lower temperatures at which leptogenesis takes place. Thus, 
the production of unwanted relics, most notoriously of gravitinos within supersymmetric 



models [llHl3|. along with the lepton asymmetry asymmetry can be avoided. A lower 



energy scale might also give rise to new experimentally accessible signals connected with 



leptogenesis, see e.g. [14H16 



However, the origin of the masses of the right-handed neutrinos is yet unknown and 
their masses may well be neither hierarchical nor degenerate. When the hierarchical 
limit Mi <C M 2 no longer applies, it is well known that the decay asymmetry of Ni 
is not simply proportional to M±/M 2 with negligible corrections, as can be verified by 
inspection of the vertex and wave-function contributions to the decay asymmetry of 
iVi [H, 0]- In the finite-temperature background, there are additional corrections due 
to new cuts. While in the vacuum background, the CP- violating contribution from 
the vertex function arises exclusively from the cut through {£,4>}, where the internal 
lepton and Higgs boson are on-shell, at finite temperature also the two other possible 
cuts through {£, N 2 } or {<(),N 2 } in the vertex diagram contribute, cf. Figured] (B). 
This is because in the finite-temperature background, the cut-particles do not need to 
correspond to stimulated (suppressed) emission processes for bosons (fermions), but they 
can also correspond to absorption processes of particles from the plasma. The presence 



of these cut contributions has been mentioned in Ref. [17(. However, by now, neither 
analytical expressions for these terms have been derived nor have these been evaluated 
numerically in order to compute effective decay asymmetries or the resulting lepton 
asymmetry. These are the main goals of the present work. 

The new cut contributions are a finite density effect, and a powerful method of de- 
scribing out-of-equilibrium field theory is given by the Schwinger-Keldysh Closed-Time- 
Path (CTP) formalism 18|, ll9(. This approach has been applied to leptogenesis and has 
resulted in some recent activity which we build upon within the present work 2CH27 . 
Main advantages of the CTP approach to leptogenesis over the conventional description 
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Figure 1: Diagram (A) represents the vertex correction to the lepton self-energy in 
the CTP formalism. Diagram (B) is the subdiagram of (A) that accounts for the decays 
and inverse decays of the out-of-equilibrium particle N%. We indicate the various cuts 
that arise from demanding that the cut particles in the loop are on shell. The solid circle 
represents the standard cut through {£, <$} that is the only contribution in the vacuum or 
when finite density effects are neglected. The dashed cut is the contribution for off-shell 
t through {<p,N 2 }, the dotted cut the the contribution for off-shell <j) through {£, N 2 }. 



by semi-classical Boltzmann equations may be seen in the absence of the need of an ex- 
plicit subtraction procedure for real intermediate states (RIS) [28[ and in the systematic 
treatment of finite-density corrections. 

Within the CTP approach, the vertex diagram in Figure Q] (B) appears as a subdia- 
gram in the self-energy Figure Q] (A), which is a contribution to the lepton self-energy. 
This self-energy in turn enters the collision term of the Kadanoff-Baym equation for the 
lepton, that can be reduced to a kinetic equation which describes the gain and the loss 
and therefore the time evolution of the lepton number density. In order to simplify the 
collision term to a manageable form, it is useful to substitute equilibrium propagators for 
I and and to employ Kubo-Martin-Schwinger (KMS) relations. The present work relies 
strongly on Ref. [5[, where these simplification strategies are explained and justified in 
detail and where also many definitions and quantities that we use here are introduced. 
An additional simplifying assumption that we introduce here, but that will not hold 
true in general, is that also is maintained in equilibrium (e.g. through interactions 
with an additional lepton flavour), such that N\ is the only out-of-equilibrium particle. 
We leave a study of the situation where more than one of the right-handed neutrinos 
deviates from equilibrium to future work. 

Provided the initial distribution function for N\ is thermal, as we assume in the 
present work, the main contributions to the lepton asymmetry occur at times when the 
value of z = Mi/T is in the range of around one up to a few, for an in-detail discussion, 



see Ref. |29|. This is because at these temperatures, Ni becomes non-relativistic and 
therefore deviate from equilibrium. As a consequence, we anticipate the new cuts to be 
important under the following conditions: 
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Since the effect is due to the finite densities, M 2 should not be much larger than 
Mi, because otherwise the new contributions are Maxwell-suppressed, just as the 
equilibrium distribution functions for energies of order M 2 much larger than T. 

The new effects will be most pronounced in case leptogenesis occurs at comparably 
low values for z, when it is around one. This is the case in the transitional regime 
between strong and weak washout. In such a situation, the finite density effects 
that include the contributions from the new cuts, take their largest relevance. 

A loophole to these arguments are situations where Mi > M 2 . (We denote within 
this work by Aq the neutrino that deviates from equilibrium, whereas N 2 is assumed 
to be very close to equilibrium at the time relevant for leptogenesis. This definition 
differs from what is commonly used in the literature, where Aq corresponds to the 



lightest right-handed neutrino, cf. Refs |3CH32l|. In this case, the finite density 



effects from the N 2 are also important when leptogenesis occurs at larger values of 
z, since the distribution of N 2 can still be unsuppressed. 

The plan of this paper is as follows: In Section [21 we give the expressions for the CP- 
violating source terms that bias the lepton number and that are valid for finite ratios of 
M\jM 2 . For the standard cut through we quote the result from Ref. 0, while 

for the new cuts through {£, N 2 } and {<fi, N 2 }, respectively, we derive new expressions. 
In Section [3], we present the first numerical results for the cut through {£, <$} at finite 
Mi/M 2 and finite density, as well as for the new cuts through {£, N 2 } and {<j), N 2 }. We 
first define and evaluate expressions for the effective CP violation for several ratios of 
Mi/M 2 . Then, we solve the kinetic equations for the lepton asymmetry for the same 
values of Mi/M 2 and specific illustrative choices of Yi, Y 2 . The results for the lepton 
asymmetry are compared with the effective CP violation as well as with what is stated 
above on the parametric regions where we anticipate the new cut contributions to be 
most relevant. Within flavoured leptogenesis, a new cut in the wave- function of the 
lepton i contributes to the asymmetry. A rough estimate of this effect is presented in 
Section HI and it is found to be negligibly small. We summarise and conclude in Section [5j 



2 Thermal Vertex Function 

2.1 Vertex Self Energy and Collision Term 

Within the present work, we extend the results of Ref. . There, the kinetic evolution 
equation for the lepton asymmetry is expressed as 

^-(n E -n e ) = W + S, (1) 
dr] 

where 77 is the conformal time in the Friedmann background, W the washout term, S the 
source term and ni (fig) the comoving (anti-)lepton number density. In the radiation- 
dominated Universe, the scale factor is given by a = a^rj, where or is an arbitrary 
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constant. The washout term W is discussed in detail in Ref. [5[. It encompasses the 
tree-level decay and inverse decay processes of the right-handed neutrino. The source 
term accounts for the CP-violating loop effects. It decomposes as 



where CJ 1 is the wave- function contribution to the collision term, as it is given in Ref. [51]. 
In this work, we are primarily concerned with the vertex contribution CJ and the new 
corrections that it acquires when compared with Ref. jij . This term can be expressed in 
the usual Kadanoff-Baym form 

— tr ^f>(k)P L iSm - ^f<(k)P L iS>(k)\ , (3) 

for X = v. The Wightman self-energy ^Z v> has the diagrammatic representation of 
Figured] (A). It is given by 0] 

{iSU~P)C [iSj(p + k + q)]' ChS T Nj (-q)iA<(-p - fc)iAj(-g - k) 
-iSU-p)C [iS<(p + k + q)f CHS T N] (-q)iA<(-p - k)iA<(-q - k) 
-iS> t (-p)C [iS>(p + k + q)f ChS^i-q^i-p - fe)iAj(-g - fc) 
+i4 l (-p)C ( [iSf (p + k + Cti5->.(-g)iAj(-p - k)iA<(-q - k) 



where we sum over the indices The self-energy ^Z v< follows when applying the 

replacements <-H-> and T -h- T. For the propagators S 1 ^, and SVi, we use the 
zero- width approximations as written down in Ref. [5|. 

To be specific, let us now consider the term with % = 1 and j '• = 2. The CP- violating 
contributions from the decays and inverse decays of Ni arise when two out of the three 
propagators iSn{p + k + q), iA^—p — k) and iS^i—q) are on shell. These are the 
cut propagators. As a consequence of this, the third propagator is off shell. Since the 
> and < propagators are purely on shell, it follows that only terms where the off-shell 
propagator is of the time-ordered T or anti-time-ordered T type contribute. The collision 
term ([3]) can therefore be split into the portions 

CJ = C^ + C^ N2 + CJ m \ (5) 

where the superscripts indicate through which two of the loop propagators in the vertex 
correction the cut goes. Likewise, we decompose the vertex contribution to the source 
term as 

Cjiv ctvicj) _|_ c<v(f>N2 _|_ (rviN2 ^g-^ 
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Within the present work, we make the simplifying assumption that N 2 is in thermal 
equilibrium at all times, whereas Ni is in equilibrium initially but then deviates from 
equilibrium when it becomes non-relativistic. We note that in general, both Ni and N 2 
will deviate from equilibrium at the same time. In fact, situations are conceivable where 
the initial abundance of N 2 is zero and Y 2 is so small that N 2 does not equilibrate before 
becoming non-relativistic. Therefore, the equilibrium deviation of N 2 can be large when 
compared to Ni, and inverse decays of N 2 may largely enhance the lepton asymmetry. 
For now, leave this interesting possibility for future work and note that N 2 may be 
maintained in equilibrium by a stronger coupling to a different lepton flavour within 
which no asymmetry is produced. 



2.2 Cut through {£, 0} 

The term C 1 /^ arises from those terms within i^J <,:> , Eq. Pjh where the propagators iS^'J 
occur. It can be further simplified when approximating the distribution functions of I 
and 4> by Fermi-Dirac and, respectively, Bose- Einstein equilibrium distributions /^ q (k) 
and /^ q (k). Accounting for the fact that £ deviates from equilibrium only leads to higher 
order corrections in the gradient expansion, as it is explained in Ref. [5|. The gradient 
expansion corresponds to an expansion in powers of Yj 2 or equivalently H/T, where H 
denotes the Hubble rate. 

The contribution of the cut through {£, 0} to the source term can be factorised into 0] 



r*f * v"- Mi Mo) - [ d ' k ' d " k " (2n)Wk - k' - k") k> MlM2 
1„ [k,p ,M U M 2 ) - J (27r)32 | k ,| (27r )3 2 |k"| (27rj 6{k k k j/ > (A;'-p") 2 -M 2 2 



x 



f) 2 - M| 

(7) 

[1 - /; q (k') + /; q (k'o] 



and 



v»(k,M u M 2 )=J (27r )32 |p1 (2 ^2 |p1 (^( k -P'~ p") r?ft P"; M U M 2 ) 

(8) 



x 



[1 - /; q (p) + n q (p" 



with the result 

s ,« = J 0_ cr(p , } = 4 lmKYf] J _ j|_ Sfm(k) ^ MhM2) 

(9) 

The term 5fNi{k) denotes the deviation of fNi(k) from the equilibrium Fermi-Dirac 
distribution, Sf N1 (k) = f N1 (k) - /^(k). 

In Ref. [5( , this result has been evaluated in the hierarchical limit Mi / M2 1 , which 
has the virtue that, as it is explained in the Introduction, the evolution of the lepton num- 
ber density is independent of M 2 up to an overall proportionality. Furthermore, when 
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M\jM 2 <C 1, the collision term (Q can be reduced analytically to a one-dimensional 
integral, which can easily be evaluated numerically. For the present work, we numer- 
ically evaluate a multi-dimensional integral for S ve< ^ that remains when exploiting the 
5-functions. We note that also within the related articles 22|, |24[, all source terms have 
been evaluated in the hierarchical limit, such that here, we present the first quantitative 
results for effects of finite M\jM 2 in a finite-density background. 

2.3 Cut through {0, N 2 } 

Now the cut goes through the internal N 2 and <p lines of Figured] (A), and the internal 
i is off-shell. The according contributions to the collision term are 



w<j>N2, 



Y *2 Y 2 f dk ° d4 P A 
i j J 2vr (2tt) 4 (2tt) 4 



(10) 



tr{iS<(k)iS> t ( 


-p)C [Sj(p + k + q) 




-?)iAJ( 


-p- 




-9 




-iS<(k)iSU 


-p)C [Sj(p + k + q) 




-9)iAj( 


-p- 


fc)iAJ( 


-9 


-*) 


+iS>(k)iS< i ( 


-p)C [Sjip + k + q) 




-g)iAJ( 


-p- 


*)iAj( 


-9 


-*) 


-iS>(k)iSU 


-p)C [Sj(p + k + q) 




-9)iA?( 


-p- 




-9 


-*)} 



As explained above, we assume here that besides I and </>, N 2 is in equilibrium, such 
that its number density is given by the Fermi-Dirac distribution function /^ 2 (k). This 
allows us to apply KMS relations and two further replacements, which become identities 
under the integrals (c/. Ref. [5|, where this is explained in more detail): 



LSjf (- P )iA T /(-p -*)-►- [LS&(-p)iAj(-p - fc) + LSfo(-p)iAj(-p - fc)] 

<f (-?)iAf * (-9 - A) i [i^,(-g)iA>(-g - A) + i^(-g)iA<(-g - k)] 
The result of these simplifications is 

c v / N2 (k) = - \y*' z y 2 2 - y; 2 y* 



11) 



■2 V *2 } I dk° d 4 p d 4 g 
2 1 1 2ir (27r) 4 (27r) 4 



(12) 



1 

-tr 
2 



iSf(k)VS^(-p)C [Sf(p + k + q)YcHA<(-p - k) 
-iS>(k)iSU-p)C [Sf(p + k + q)f CHA>(-p - k) 
x \iS< 2 (-q)iA>(-q - k) - iS> 2 (-g)iA<(-g - k)] } . 
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The final contribution to the source term, 



•v<j>N2 



I 



d 3 p' 
(2^)3 



(13) 



then follows when substituting the explicit forms of the propagators (as they can be 
found in Ref. I5|) as 



To our knowledge, this is the first report of a result for a source term for leptogenesis that 
is present at finite density but absent in the limit of a vacuum background. We remark 



that the last factor results from the expression —f^ — fm = ~ [^- Jr fl q ]fm~ fd> q [l~ fm\- 



This could also be guessed starting from a hypothetical loop factor when the decay 



I -»■ 4>N 2 was kinematically allowed, [1 + f cq - f$ 2 ] = [1 + /^ q ][l - f^ 2 ] + fl q f^ 2 , and 



applying the replacements — f^ 2 -H- [1 — f^ 2 ]- This argument may be considered as a 
consistency check for our derived result (I14p . 

It is important to notice that for M 2 ^> T, the source term (|T4l) is strongly Maxwell 
suppressed. While in such a situation, /^ 2 (k') is always suppressed because of the large 
mass of N 2 , the energy- momentum conserving (^-functions always imply that then at 
least one of the momenta p or k" is much larger than T, such that f^ q (k") or Sf ni(p) 
are suppressed as well. 



We finally consider the cuts through the propagators N 2 and t within the loop and take 
(f) to be off-shell. The contribution to the collision term is 




2.4 Cut through {£, N 2 } 




(15) 



p)C [S>(P + k + q)f C*iSfo(-g)iAj(-p - fc)iAj(-g - k) 
p)C [Sj(p + k + q)] t ChS^i-q^i-p - k)iA< (-q - k) 
■p)C [S<(p + k + q)f C*LSfo(-g)iAj(-j, - fc)iAj(-g - k) 
p)C [Sj{p + k + q)]' ChS^(-q)iAl(-p - k)iA>{-q - k)} 
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Again, we substitute for Sn2, Sg and the equilibrium propagators and employ KMS 
relations. Under the integrals, we make the replacements 

i5 jvf (-p)iC S e' T (P + k + q) 



C f (16) 

iSU~P)iC [Sf(p + k + q)f + iSU~P)iC [S<(p + k + q)f Ct" . 
Furthermore, it is useful to notice that when substituting the equilibrium propagator for 



1 

^2 



LS^AOLS&C-iOC [5>(p + fc + g)f C^-^) 
-iSf (A0lS&(-p)C [S<(p + fc + q)f ChS^(-q) 

is odd under the exchange k°,p°, q° — > —k°, —p°, —q°, while Im[iA^' T (— q — k)] (which is 
the off-shell contribution) is even. Making use of these additional remarks, the collision 
term simplifies to 



C 



wlN2 



[Y *2y2 _ y2y^ J 



dk° d 4 p d 4 q 



2tt (2vr) 4 (2?r) 4 
tr{ [LSfo,(- g )LST(AO - iS> 2 {-q)iSf{k)] iAj(-p - k) 
x [i5^(-p)C [52= (p + * + q)f CHA<(-q - k) 
-\S> Nl {-p)C [S>(P + k + q)f CHA>(-q - k)] } . 



(17) 



The source term is 



,5' 



v£N2 



d 3 p' rv eN2( ri > 
rL'A \ P 



(2tt) 3 * 

which becomes, when substituting the finite-density propagators, 

d 3 k" d 3 p d 3 p" d 3 k' 



S viN2 =Alm[Y 1 %* 2 ] J 



(18) 



d 3 p' 



(2vr)32|k"| (2^)32^ + M 1 2 (2tt)32|p"| ( 27r ) 3 2 ^k' 2 + M| (2vr) 3 2|p'| 

^M 1 M 2 (19) 



x (2tt) 4 <5 4 (p - fc" - /)(2tt) 4 5 4 (A; / - fc" - p') ^ 



(p" + k') 

xsfm(p) [i - /rv) + /; q ( k ")]] >< t/rv) - /&ooi . 

Again, we note that the factor [fP — f^ 2 ] = — [1 — /^ q ]/^ 2 + — fm] could also be 
guessed from a would-be factor [1 — /^ cq ][l — f^ 2 ] — fPf^ = 1 ~ fPfmi which would 
arise if — > £N 2 was kinematically allowed, and the replacements — f^ 2 -H- [1 — f^ 2 ]- 
We also note that S veN2 is Maxwell suppressed for M 2 3> T, in analogy with what is 
discussed S V ^ N \ 
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3 Examples 



We now present results from the numerical evaluation of the source terms (I9jl4lli9p . In 
Section 13.11 we define an effective measure for the CP violation from the various cut 
contributions through a benchmark out-of-equilibrium distribution for Nx. Then, we 
evaluate this effective CP violation as a function of z — M\/T for various values of 
Mi/M 2 . We find that in case M 2 > M\ the new corrections are only significant for z ~ 1 
or smaller. If M 2 < M±, the corrections from the new cuts can be relevant for larger 
values of z. We proceed in Section 13.21 with the calculation of the lepton asymmetry 
by solving the kinetic Boltzmann-type equations ([T]) with the new source terms. Our 
choice of the washout strength is motivated by the wish to exhibit models where the 
contributions from the new cuts are sizable, that is where relevant contributions to the 
final asymmetry arise at values for z around one. This is the case in the transitional 
regime from weak to strong washout. We find that indeed, the new cut contributions 
can have a sizable impact on the final asymmetry. However, in case sizable contributions 
to the asymmetry originate from z <C 1 in our simulations, the quantitative results need 
to be interpreted with care. This is because in these regions, thermal corrections to the 
masses and widths of </> and i will become relevant. 

3.1 Effective CP-Violating Parameter 

In order to obtain a quantitative comparison of the amount of CP violation from the 
various source terms, it is useful to define a benchmark form for the distribution <5/jvi(k). 
This is necessary, since the precise form of 5/jvi(k) depends on time, washout strength 
and initial conditions. We follow Ref . j^] by taking for f^i (k) a Fermi-Dirac distribution 
with a pseudo-chemical potential /ijvi- The deviation from equilibrium is then obtained 
by expanding to linear order in fi^i/T, 



When the only interactions of the neutrino Ni are mediated by Y\, the actual distribu- 
tion function is not exactly described by the pseudo-chemical potential, and Eq. ([20]) 
should indeed only be regarded as a useful benchmark for the purpose of comparing the 
various contributions to the source term. Note however that in case there are fast elastic 
scatterings between the neutrinos N±, Eq. ( 1201) should be a very accurate description 
for the actual distribution. We substitute Eq. ([20} into Eqs. f[9lfTlfT9]) and into S wi and 
Sm2>mi as given in Ref. (if, where Sm2>mi is the source term including both, vertex 
and wave function terms, evaluated in the hierarchical limit. Out of these, we take the 
ratios 



*/™(k)=/£(k)(l-/£(k)) 



A* ATI 



(20) 



(21) 



and 



(S v + S wi )/ SW 2 >Afi , 



(22) 
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as functions of z — M\/T and with S v defined in Eq. (Q. The ratio (I2ip allows for a 
comparison of the source term at finite density and the standard cut through £ and <fi 
only with the source term in the hierarchical limit. Through the ratio f j22l) . we compare 
the source with all cuts at finite density with the hierarchical limit. To quantify the 
effect of the new cuts, both ratios (|21[) and (]22[) are compared with one another. 

Before we present the numerical results for Eqs. ( T2"Tj) and (122]) . a few remarks on their 
relevance and range of validity are in order. First, note that the relevant contributions to 
leptogenesis are generated when z ~ 1 or larger. (This is to be understood in the sense 
that while z is much smaller than one, no sizable contributions to the final asymmetry 
are generated). This is a consequence of the fact that N\ must become non-relativistic 
and Maxwell suppressed before it equilibrates, because otherwise the lepton asymmetry 
would be completley washed out. Therefore, the effect of the new contributions is more 
relevant if it extends to larger values of z. 

Second, we remark that that for small z, finite temperature effects should become 
relevant, which we do not take into account within the present work. The most important 
correction is the contribution of the top quarks to the Debye mass-square of the Higgs 
boson, m 2 H = (l/4)/ijT 2 , where ht is the top-quark Yukawa coupling. Besides, also the 
SU(2)l gauge couplings of the Higgs bosons and the leptons are of relevance. A full 
evaluation of these finite-temperature effects in the context of leptogenesis has not been 
performed yet, but it would be of great importance for obtaining quantitatively accurate 
results in models where asymmetries generated at values of z that are somewhat smaller 
than one are relevant. Moreover, the size of the top-quark Yukawa coupling at the 
energy scale of leptogenesis depends due to running on the mass of the Higgs boson, 
which is yet unknown. For the present discussion, we therefore keep in mind that for the 
contributions that are generated at z < 1, we have to anticipate an inaccuracy of order 
one. 

In Figure El the ratios ( 12T|) and ( 122]) are plotted for several values of Mi/M 2 . The 
main features for each of the particular values can be summarised as follows: 

(A) : For M 2 /Mi = 1.1, the hierarchical limit is of course not a good approximation. 

The total asymmetry is dominated by the resonantly enhanced contribution S . 
As a consequence of this, (S ve< ^ + S wf )/ Sm 2 ->m 1 is much larger than one. Since the 
distribution of N 2 is only very weakly thermally suppressed when compared to Ni, 
gvicj) _|_ gwf an( j gv _|_ gwf a g ree on iy f or comparably large values of z. 

(B) : For M 2 jM\ = 2.0, we note that the hierarchical limit becomes a better approxi- 

mation to S ^ + S wi , while yet deviating by about 20%. Compared to case (A), 
S ^ + S wf and S v + S wi start to agree for smaller values of z, as the distribution 
of N 2 now suffers from stronger thermal suppression. We note that the very sub- 
stantial deviation of S vi(p + S wi from S v + S wi for very small values of z needs to 
be interpreted with care in the light of the thermal corrections mentioned above. 

(C) : For M 2 jM\ = 5.0, the results obtained from the hierarchical limit Sm2>Mi and 

the source S ^ + S wf with the standard cuts only are now in good agreement, as 
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(A): M 2 /M 1 = 1.1 



(B): M 2 /M 1 = 2.0 




Figure 2: {S v + S wi )/S M2M (solid line) and (S v ^ + S wt )/S Ma>Ml (dashed line) over 
z = Mi/T. 

anticipated. This serves also as a consistency check for our numerical evaluation 
of S v + S wi in the form of a multi-dimensional integral as compared to Sju 2 >A/i 
through the one- dimensional integral given in Ref. [5[. The new cut contributions 
become suppressed for even smaller values of z along with the stronger Maxwell 
suppression of N 2 . The comment regarding Scenario (B), that the deviations that 
originate for values z <C 1 need to be interpreted with care due to thermal correc- 
tions, is even more relevant for the present case. 

(D): In this Scenario, M 2 jM\ = 0.5. Since N 2 is now lighter than JVi, the coupling Y 2 
must be chosen sufficiently small, such that N 2 does not equilibrate and wash out 
the lepton asymmetry before it becomes non-relativistic and Maxwell suppressed. 
In this situation, for the whole relevant range of z, neither the hierarchical limit 
nor the sources S ^ + S wi are an accurate approximation to the full result S v + S . 
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3.2 Lepton Asymmetries 



From the comparison of the effective amount of CP violation in presence of the new 
cuts with the case when the new cuts are neglected, we see that when M 2 > Mi, the 
deviations are only sizable for z is of order a few or smaller. For larger values of z, the 
distribution of the is strongly Maxwell suppressed, and the new cut contributions 
become irrelevant. In order to exhibit the effect of the new cuts, we therefore choose 
the parameters Mi and Y± such that leptogenesis takes place in the transitional regime 
from weak to strong washout. In this situation, sizable contributions to the final lepton 
asymmetry are generated for z ~ 1, such that the effects of the new cuts becomes 
relevant. 

We choose for N\ )2 thermal initial conditions. In the early Universe, these may 
be established through interactions via heavy gauge bosons that freeze out at times 
before leptogenesis takes place. We assume that N 2 is maintained in equilibrium due 
to a Yukawa coupling with an additional lepton flavour, within which no asymmetry is 
produced. The coupling Y 2 is chosen smaller than Y 1; such that washout effects from 
inverse decays of N 2 are negligible. (Explicitly, for the scenarios with M 2 > Mi, the 
largest error from washout through Y 2 occurs for M 2 jM\ = 1.1. Since the thermal 
suppression of N 2 compared to Ni is very small in this case and we have chosen Y 2 /Yi = 
1/2, we expect that the washout is underestimated by 20%. The accuracy improves for 
larger ratios of M 2 > Mi due to the thermal suppression of N 2 and its irrelevance for 
washout). While this appears as a somewhat special setup, we note that even when Y 2 
is large, the vector of couplings of N 2 to the various left-handed lepton flavours defines 
a particular linear combination £ 2 of leptons that are washed out through inverse decays 
of N 2 . This linear combination can in general be linearly independent of the linear 
combination I within which the lepton asymmetry through decays and inverse decays 
of Ni is produced. The contribution to the asymmetry in £ that is orthogonal to £ 2 is 
then unaffected by the washout through N 2 [3l|. (Note the different assignment of the 
heavy neutrinos to the indices 1,2 in that work). Therefore, the qualitative features 
of our particular setup should be relevant for parametrically more generic models of 
leptogenesis. However, it would still be interesting to study the effect of different initial 
conditions and the possibility of N 2 being out-of-equilibrium within future work. 

As it is described in detail in Ref. [5[, we obtain the numerical results as follows: 
First, we solve the evolution equations for /Wi^k) as a function of z. These, we feed into 
the washout term W and the source term S for the leptons, in order to solve for the 
lepton asymmetry in Eq. ([T]). The expansion of the Universe is taken into account when 
inserting the scale factor according to Mi 2 — > aM li2 . 

The results for the lepton-number to entropy ratio Yg as defined in Ref. j^] are pre- 
sented in Figure [3j Note that it is instructive to compare the particular panels with 
those in Figure |2j Again, we summarise some features for each of the particular values 
of M 2 /Mi: 

(A): For M 2 /M\ = 1.1, the hierarchical limit clearly underestimates the full result. 
This is anticipated, because the wave-function contribution to the full result is 
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Figure 3: Evolution of the lepton asymmetry Yi over z = M\jT . The choice of parame- 
ters is Mi = 10 13 GeV, Y\ = 2i x 10~ 2 . Solid: full result; dashed: result with contribution 
from cut through {£,<p} only; dot dashed: hierarchical limit M 2 ^> M\. 



resonantly enhanced for Mi fa Mi- The new cuts give rise to relative corrections at 
the 10% level. Note however that the absolute correction is larger when compared 
to scenarios (B) and (C). The relative correction is marginalised due to the resonant 
enhancement of S . 

(B) : For M 2 /Mi = 2.0, there are sizable deviations of the full result from the result 

with the cut through {£, 0} only, that arise in the region z ~ 1. 

(C) : For M 2 /Mi = 5.0, the new contributions are of importance for smaller values of z 

when compared to scenarios (A) and (B). In the light of the thermal corrections 
that we anticipate to be important for small values of z, the quantitative result 
needs to be interpreted with care. The good agreement between the result from 
the standard cut through {£, <fr} only and the hierarchical limit serves again as a 
consistency check for our numerical evaluations. 

(D) : For M 2 /Mi = 0.5, we choose a smaller value for Y 2 , motivated by the requirement 

that the A^ 2 must not wash out the lepton asymmetry. The full result and the 
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Figure 4: Wave function correction that contributes to the lepton asymmetry in 
flavoured leptogenesis. 

result with the cut through {£ , 0} only receive different contributions even for 
values of z of the order of a few, cf. Figure [2] (D). We expect therefore only a small 
contamination from theoretical uncertainties due to thermal corrections. Note that 
in this scenario most of the lepton asymmetry of the Universe originates from the 
new cuts. 

4 Flavoured Leptogenesis 

When we insert a loop of <p and as a wave-function correction into the propagator 
Si in the vacuum, no branch cut term due to on-shell <fi and N 2 arises for kinematic 
reasons. Again, this holds no longer true in a finite temperature background. Therefore, 
CP-violating source terms can arise from the diagram in Figure HI Since there is no 
lepton number violation, this does not yield a contribution to the asymmetry in models 
of unflavoured leptogenesis. However, this diagram may appear as a source in flavoured 
scenarios. In this Section, we give a rough estimate of this contribution, leading to the 
conclusion that it is generically negligible. 

Within the CTP-formalism, the form of the Wightman type self-energy, that is given 
diagrammatically in Figure HI reads 




(23) 



{iS> t (p)iA<(p - k)iSj cd (k)iS T Nj (q)iAl(q - k) 
-iSfo(p)iAJ(p - k)\SUk)\S> N] (q)\A<(q - k) 
-XSUP)^(P ~ k)iS> d (k)iS T Nj (q)iAl(q - k) 
+LSj i (p)iAj(p - k)iSf cd (k)iS>(q)iA<(q - k) 



} 
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Compared to the unflavoured scenario, we have promoted the Yukawa couplings of the 
right-handed neutrinos iVj to a matrix Y ia , where the first index refers to the right-handed 
neutrino and the second index to the left-handed lepton flavour. For the definitions of 
the model Lagrangian and the lepton propagator Se, we refer to Ref. 
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Next, we insert ^™j, £> into the collision term 



. This again simplifies when sub- 
stituting equilibrium propagators for N2, I and and exploiting KMS relations. We 
furthermore choose to work in the flavour basis where the matrix of Standard Model 
lepton Yukawa couplings h a b is diagonal. This is advantageous since off-diagonal compo- 
nents of Sfab are damped away quickly in this basis, provided the interactions mediated 



by h a b are f as t compared to the Hubble rate [14j, |33H35| . cf. Ref. [26j for a description 



within the CTP approach and for a numerical study of this effect. Note that this also im- 
plies that the off-diagonal components of the equilibrium propagator for t are vanishing. 
We eventually obtain for the collision term 



L iaa\ K ) — - [ Y la r lcl2c r 2a ~ 1 2a 1 2c* l c * la\ 



dk° d 4 p d A q 



2tt (2tt) 4 (2tt) 4 

itr{ (i5&(A0LS$,(p)iAj(p - *) - tf^(*)LS&(p)iAJ (p - h)) iRe 
x (iS<^ q )iA>( q - k) - LS^faJiAJfe - k)) 



(24) 



k 2 + ie 



where e is infinitesimal. In the zero- width limit, this term contains a divergence from the 
factor l/{k 2 + ie) that originates from the propagator Sf(k), since Sf aa {k) is proportional 
to S(k 2 ). At finite temperature, this is regulated through the replacements 



1 



k 2 + i 



1 
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k 2 



mj c + ik°T £ 



and 



S(k 2 ) -+ — 



2n 



1 



k 2 



+ ik°Ti 



m 



ifcor. 



(25) 



(26) 



Here, the thermal masses of the leptons are m 2 a = h 2 aa ^ ,h {k) + g 2 <; hl,9 (k), where g is 
the SU(2)l gauge coupling. The functions ^^{k) and ? bl ' 9 (fc) are of order T 2 for k of 
order T. They account for thermal mass corrections from flavour- sensitive and flavour 
blind interactions, and are discussed in more detail in Ref. (ifj. Likewise T £ = g 2 g hl (k) 
is the finite-temperature width of the leptons, where g u (k) is of order T when k of 
order T. When leptogenesis occurs at temperatures of roughly below 10 11 to 10 12 GeV 
but above 10 8 to 10 9 GeV, the r-lepton Yukawa coupling h T is in equilibrium, while the 
electron e and muon \x couplings are yet out-of-equilibrium. In this situation, flavoured 
leptogenesis distinguishes effectively between two flavours, where flavour 1 is identified 
with t and the coupling hn = h T and flavour 2 with a linear combination of e and fi 
and with negligible Yukawa coupling h 2 2 ~ 0. The integral <HM could possibly again be 
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evaluated numerically. However, we can estimate that from the factors f )25|) and ( 126)) . 
the /c°-integration yields a factor 

m "~ m ^ „ JLS (27) 

for |k| ~ T. The latter estimate follows from above assumption that only h T appears as a 
relevant Standard Model Yukawa coupling. Above suppression factor is to be compared 
to 1/T 3 for a lepto genesis scenario with Mi ~ M 2 but no pronounced resonant enhance- 
ment. Since h^/g A <C 1, we conclude that this contribution to flavoured leptogenesis is 
suppressed due to the large width of the leptons I at finite temperature. In other words, 
the contribution to the CP- asymmetry is rendered ineffective because the separation 
between the resonances of the lepton quasi-particles is well within their overlap due to 
the finite widths. 



5 Conclusions 

In this paper, we have presented in the form of Eqs. (1141) and (1191) the first results for 
source terms that contribute to the lepton asymmetry in a finite-density background, but 
that are absent in the vacuum. In order for these to be relevant, we have seen that M 2 
should not be much larger than T at the time of leptogenesis. The main features of the 
numerical evaluations are easily understood. First, M 2 should not be much larger than 
Mi for the new effects to be sizable at larger values of z, cf. Figure [3 We reemphasise 
that the quantitative results for the scenario in Figure [3] (C) should be considered with 
great care since a substantial amount of the deviations incurred by the new effects is 
generated at very small values of z, where we expect thermal corrections to become 
relevant. Second, the effect is most pronounced in scenarios where a sizable amount of 
the asymmetry is produced for comparably small values of z, as it is the case for the 
transitional regime between weak and strong washout. As a loophole, we find that the 
largest effects arise within a somewhat unconventional scenario with M 2 < Mi, since in 
this situation the new cuts are important throughout the time of the out-of-equilibrium 
decays of Ni. Therefore, relevant contributions from the new cuts also result for larger 



values of z. However, in regard of Ref. [3l|, where it is shown that decay asymmetries 



from the heavier singlet neutrinos generically survive subsequent washout, such a scenario 
may not only appear as a mere loophole. Provided the reheat temperature is large enough 
to produce the heavier singlet neutrino in the early Universe, the present work implies 
that a large contribution to the baryon asymmetry of the Universe generically results 
from the new cuts. 



We mention that while this work has been in preparation, Ref. [36| appeared, where 
non-standard cuts are discussed as well. Since finite density effects are neglected, it is 
found there that the new cuts are only allowed within certain scattering processes. More 
precisely, in order to satisfy the kinematic thresholds for the presence of cut contributions 
in the vacuum, diagrams with additional radiation of Standard Model particles are con- 
sidered. Note that if the resulting amplitudes were substituted as source terms into the 



1(3 



kinetic equations for the lepton asymmetry, they would be subject to the same Maxwell 
suppression that occurs also in the present case for M 2 ^> T. However, compared to 
the new source terms that we have derived in the present work, the terms discussed in 



Ref. [36| are subject to further suppression because of the insertion of additional coupling 
constants. 

On the ground of the present results, it would be interesting to address the following 
points in the future: 

• The results ( fl4|) and (fl9|) should be generalised to include effects from deviations 
of the distribution f^2 from equilibrium, which we expect to be non-zero in generic 
scenarios of leptogenesis. 

• Systematic investigations of the parameter space for scenarios where the new cuts 
are relevant would be desirable. We expect the new cut contributions to be of 
crucial importance for phenomeno logical studies where it is assumed that M 2 < M±, 



such as Refs. |30H32| . where N 2 may or may not be in equilibrium. (Note the 
different definitions of in these works). Note that in case M 2 <C Mi, thermal 
corrections may become important again, since the thermal masses of (ft and i at 
the time of the decay of N\ may exceed M 2 . We expect the new cuts still to be 
important in such a situation, but the interpretations of the particular cut particles 
in terms of absorption processes may change to emissions and vice versa. 

• In the present work, we have restricted ourselves to compute the lepton asymme- 
tries for thermal initial conditions for Ni >2 - This is in part motivated by the fact 
that for a vanishing initial density of N\, the final asymmetry is a remainder of 
an incomplete cancellation of a contribution that is created initially at small z, 
when Ni is underabundant, and an opposite one through later decays when z is 
larger and N\ overabundant. The fact that the cancellation is incomplete is be- 
cause the opposite asymmetries are affected differently by the washout, since they 
occur at different temperatures. In order to obtain a quantitatively reliable result 
for the remaining asymmetry, a rather accurate prediction of the asymmetry that 
is created through inverse decays at small z is necessary. Due to the uncertainties 
because of thermal corrections, that we have emphasised in this paper, such an 
accurate prediction is presently not available. The situation somewhat improves 
for the thermal initial conditions that we consider within the present work, because 
Ni is always overabundant and the produced asymmetry is therefore of the same 
sign for all values of z. 

In regard of these points, we briefly comment on possible technical improvements that 
may prove very useful in order to increase the accuracy of the predictions for leptogenesis 
from the new cuts as well as from the standard cuts. It would be particularly interesting, 
if the following issues were addressed: 

• As it should be clear from the discussions within this paper, the uncertainties due 
to thermal corrections for small z are problematic for the predictivity of the new 
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cut contributions as well as more generically for leptogenesis in the weak washout 
regime, in particular for vanishing initial conditions for N±. In order to resolve this 
issue, a calculation of the rates Ni -H- £<f) for temperatures that are of order Mi 
or larger is necessary. Both, the tree-level rates as well as the CP-violating loop 
effects need to be calculated. For that purpose, in particular the thermal masses 
and finite widths of £ and d> should be taken into account. First work into this 



direction has been reported in Ref. [37 



• In Ref. [25], sizable effects from the finite width of N\ in the CP- violating source 
term have been reported. It needs to be verified, whether the initial conditions 
chosen in that work are applicable to the situation in the early Universe. Further- 
more, the finite widths of I and 0, which are much larger than for N\ and have 
been neglected so far, need to be taken into account in a future calculation. 

Considering these comments and a number of related papers, the present work may 
be regarded as a contribution to present efforts to improve the theoretical description 
of leptogensis, to increase the accuracy of quantitative predictions and to extend their 
range of applicability. 
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